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We present a global re-analysis of the most recent experimental data on azimuthal asymmetries in 
semi-inclusive deep inelastic scattering, from the HERMES and COMPASS Collaborations, and in 
e + e _ — V hi hi X processes, from the Belle Collaboration. The transversity and the Collins functions 
are extracted simultaneously, in the framework of a revised analysis in which a new parameterisation 
of the Collins functions is also tested. 

PACS numbers: 13.88. +e, 13.60.-r, 13. 66. Be, 13.85.Ni 

I. INTRODUCTION AND FORMALISM 

The spin structure of the nucleon, in its partonic collinear configuration, is fully described, at leading- twist, 
by three independent Parton Distribution Functions (PDFs): the unpolarised PDF, the helicity distribution 
and the transversity distribution. While the unpolarised PDF and the helicity distribution, which have been 
studied for decades, are by now very well or reasonably well known, much less information is available on the 
latter, which has been studied only recently. The reason is that, due to its chiral-odd nature, a transversity 
distribution can only be accessed in processes where it couples to another chiral-odd quantity. 

The chiral-odd partner of the transversity distribution could be a fragmentation function, like the Collins func- 
tion p] or the di-hadron fragmentation function [2-4] or another parton distribution, like the Boer-Mulders [5] 
or the transversity distribution itself. A chiral-odd partonic distribution couples to a chiral-odd fragmenta- 
tion function in Semi-Inclusive Deep Inelastic Scattering processes (SIDIS, £N — > £hX). The coupling of two 
chiral-odd partonic distributions could occur in Drell-Yan processes (D-Y, pN — > £ + £~X) but, so far, no data 
on polarised D-Y is available. Information on the convolution of two chiral-odd fragmentation functions (FFs) 
can be obtained from e + e~ — > h\ I12 X processes. 

The u and d quark transversity distributions, together with the Collins fragmentation functions, have been 
extracted for the first time in Refs. [6, 7 , from a combined analysis of SIDIS and e + e _ data. Similar results on 
the transversity distributions, coupled to the di-hadron, rather than the Collins, fragmentation function, have 
been obtained recently [5j . These independent results establish with certainty the role played by the transversity 
distributions in SIDIS azimuthal asymmetries. 

Since the first papers [3 [7], new data have become available: from the COMPASS experiment operating on a 
transversely polarised proton (NH 3 target) [HI [TO]; from a final analysis of the HERMES Collaboration [TT] and 
from corrected results of the Belle Collaboration 12 . This fresh information motivates a new global analysis 
for the simultaneous extraction of the transversity distributions and the Collins functions. 

This is performed using techniques similar to those implemented in Refs. Oil]; in addition, a second, different 
parameterisation of the Collins function will be tested, in order to assess the influence of a particular functional 
form on our results. 

Let us briefly recall the strategy followed and the formalism adopted in extracting the transversity and Collins 
distribution functions from independent SIDIS and e + e~ data. 

A. SIDIS 

We consider, at 0(k±/Q), the SIDIS process tp^ — > £' h X and the single spin asymmetry, 
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4 sm(0 h +0 s ) _ „ J d<j> h d(j> s [d<r T - do 4 ] sin(<fa + g ) 

UT " Jd4> h d<j> s [d*t + d*i] 

where da 44 is a shorthand notation for 

da™ 



g?ie e?y dz d 2 PT d(f>s 
and x, y, z are the usual SIDIS variables: 

Q 2 _ (P ■ q) _ Q 2 _ _ {P-Ph) m 

We adopt here the same notations and kinematical variables as defined in Refs. [T3], to which we refer for 
further details, in particular for the definition of the azimuthal angles which appear above and in the following 
equations. 

By considering the sm(<^ + <f>s ) moment of Aut [14] . we are able to single out the effect originating from 
the spin dependent part of the fragmentation function of a transversely polarised quark, embedded in the 
Collins function, A N D h / q t(z,p±) — (2p±/ z m^) H^~ q (z,p±) [15], coupled to the TMD transversity distribution 
A r r/(.r./,-. ) 0: 



^2 e lJ d 4>h d<p s d 2 k ± A T q(x, k±) A N D h/qf (z,p±) sin(0 s + ip + <j%) sm{(j> h + <j> s ) 

,sin(<p h +4> s ) _ _9 fo\ 

UT ~ v-^ f da > W 

22 e l J d(j} h d4>sd 2 k±_ f q /p(x,k±) — D h/q (z,p±) 



where p ± = Pt — zkj_, and 

da _ 2ira 2 2 d{Aa) _ da 1 ^ ^ da 1 ^^ _ Ana 2 

dy sxy 2 dy dy dy sxy 2 

The usual integrated transversity distribution is given, according to some common notations, by: 

A T q(x) = h u (x) = Jd 2 k 1 _ A T q(x, k X ) ■ (5) 

This analysis, performed at 0(k±/Q), can be further simplified adopting a Gaussian and factorized parame- 
terization of the TMDs. In particular for the unpolarized parton distribution (TMD-PDFs) and fragmentation 
(TMD-FFs) functions we use: 

e -k 2 ± /(k 2 ± ) 

f q / P (x,k ± ) = f q / P {x)—^y- (6) 



e -pi/{pi) 



Dh/q(Z,P±) = P>h/q{z) _ lin2 x , (7) 



with (kj_) and (pj_) fixed to the values found in Ref. [T5J by analyzing unpolarized SIDIS azimuthal dependent 
data: 

(fci) = 0.25 GeV 2 (pi) = 0.20 GeV 2 . (8) 

The integrated parton distribution and fragmentation functions, f q / p (x) and D h / q (z), are available in the 
literature; in particular, we use the GRV98LO PDF set [T7] and the DSS fragmentation function set [IS] , 

For the transversity distribution, ATq(x, k±), and the Collins FF, A N D h / q ^(z,p±), we adopt the following 
parameterizations [6j: 

A T q(x,k ± ) = - (x) [f q/p (x) + A g (a;)] 6 (9) 
A N D h/q ,(z,p x ) = 2Mq (z) D h/q (z) h(p ± ) 6 , (10) 
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with 

Nq (x) = m ^ Q (i - *) p [ a Z P (11) 



( Z ) = N? z<{\ - z) s n+ J gS (12) 
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h(p ± ) = V2e^e-rt/ M Z, (13) 

and —1 < N% < 1, — 1 < iV^ < 1. We assume (k 2 _) T = (k±). The combination [f q / p (x) + Aq(x)], where Aq(x) 
is the helicity distribution, is evolved in Q 2 according to Ref. [15]. Notice that with these choices both the 
transversity and the Collins function automatically obey their proper positivity bounds. A different functional 



form of Afq (z) will be explored in Section II B 



Using these parameterizations we obtain the following expression for 

W h j£ ^ Wf e ~7p^ )0 E 4 W (») [/«/*(*) + W (*) D h/q (z) 

A7 T iMs) = d2//d2 q ^ , ' ■ -i- 1 



-Pt/(p1) 



[i + (i-y) 2 ] 



(P|) sxy 



/?/},(«) D h/q {z) 



with 



When data or phenomenological information at different Q 2 values are considered, we take into account, at 
leading order (LO), the QCD evolution of the integrated transversity distribution. For the Collins FF, A Du 9 t, 
as its scale dependence is unknown, we tentatively assume the same Q 2 evolution as for the unpolarized FF, 

D h/q(z)- 



B. e e — > h\h,2 X processes 



Remarkably, independent information on the Collins functions can be obtained in unpolarized e + e~ processes, 
by looking at the azimuthal correlations of hadrons produced in opposite jets [2U]. This has been performed by 
the Belle Collaboration, which have measured azimuthal hadron-hadron correlations for inclusive charged pion 
production, e + e~ —> tttt X [T2J [ST] [55] . This correlation can be interpreted as a direct measure of the Collins 
effect, involving the convolution of two Collins functions. 

Two methods have been adopted in the experimental analysis of the Belle data. These can be schematically 
described as (for further details and definitions see Refs. [6] [20] , [22 ] ) : 

i) the "cos(( / 9i + ipi) method" in the Collins-Soper frame where the jet thrust axis is used as the z direction and 
the e + e~ — > qq scattering defines the xz plane; (pi and (p 2 are the azimuthal angles of the two hadrons around 
the thrust axis; 

ii) the "cos(2<^o) method", using the Gottfried-Jackson frame where one of the produced hadrons (h 2 ) identifies 
the i direction and the xz plane is determined by the lepton and the hi directions. There will then be another 
relevant plane, determined by z and the direction of the other observed hadron hi, at an angle ipo with respect 
to the xz plane. 

In both cases one integrates over the magnitude of the intrinsic transverse momenta of the hadrons with respect 
to the fragmenting quarks. For the cos(<pi + (p 2 ) method the cross section for the process e + e~ — > hi hi X reads: 

( £ a .e + e~-¥h 1 h 2 X 

dz\ dzi d cos 9 d(ipi + ipi) 

= 17 J2 e l{( 1 + cos2e ) D 'n/M) D h 2 /M) 

+ s J^l cos (<p 1+l p 2 ) A N D hl/qt {zi) A N D h2/qf {z 2 )} , (16) 
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where 9 is the angle between the lepton direction and the thrust axis and 

A N D h/q ,(z) = J d 2 Px A N D h/q t(z,p x ). 



(17) 



Integrating over the covered values of and normalizing to the corresponding azimuthal averaged unpolarized 
cross section one has: 



Rl2{zi,Z 2 ,(fil + if 2 ) 
1 



1 



1 



Sill 



(da) dzi dz 2 d((pi + 1P2) 

J2 q e 2 q A N D hl/qf ( Zl )A N D h2/qf (z 2 ) 



4(1 + cos 2 1 
2 a\ 



C0S(</31 + if 2 ) 



T, q e% D hl/q ( Zl )D h2/q (z 2 ) 



1 + 



1 



(sm 



4(1 + cos 2 



cos(<£>i + (p 2 ) P(z ll z 2 ) 



For the cos(2</?o) method, with the Gaussian ansatz (10), the analogue of Eq. (18) reads 



Ro(zi,z 2 ,tp ) = 



1 da e 



¥h 1 h 2 X 



(da) dz\ dz 2 dipo 



= 1 



1 Z\ z 2 ( 



sm 2 9 2 ) 



1 + - 



TT z{ + zt, (1 + cos 2 6 
1 z x z 2 (sin 2 6 2 ) 



cos(2(/3 ) 



J2 q elA N D hl/q ,( Zl )A N D h2/q ,(z 2 ) 



ir z\ + z\ (1 + cos 2 6 2 ) 



T, q e 2 q D hl/q ( Zl ) D h2/q (z 2 ) 
cos(2ip ) P(z 1 ,z 2 ) , 



where 62 is now the angle between the lepton and the /12 hadron directions. 
In both cases, Eqs. (18) and (19), the value of 

/cir.2 a\ 

= C(6) 



(1 + COS 2 I 



(18) 



(19) 



(20) 



can be found in the experimental data (see Tables IV and V of Ref. [2"2]). 

To eliminate false asymmetries, the Belle Collaboration consider the ratio of unlike-sign (t: + tt~ + tt~tt + ) 
to like-sign (tt + tt + + ir~n~) or charged (tt + tt + + 7t + it~ + tt~it + + n~n~) pion pair production, denoted 
respectively with indices U, L and C . For example, in the case of unlike- to like-pair production, one has 



and 



R L 12 



rt, 



1 + -C(9) cos(^i +tp 2 )Pu 
l+ l -C(6) cos((^ 1 +^ 2 )P L 



1 z\ z 2 



C(9) cos(2^ ) Pu 



1 zi z 2 



C(6) cos(2(^ ) Pl 



l + -C(6) cos( Vl + <p 2 ) (Pu - P L ) 



UL 



1 + cos(^i + tp 2 ) AH 



1 + 1 C(9) cos(2^ ) (P v - P L ) 

TT Z f + Z 2 



1 + cos(2<£>o) ^0 



and similarly for Ri 2 /Ri 2 and /Rq ■ Explicitely, one has 



Pr, 



E,e 2 [A"D n+/qt ( Zl ) A N D n - /qt (z 2 ) + A N D n - /q ,( Zl ) A N D n+/qf (z 2 )} _ ( P} 



UJN 



e l [P>K+/ q (zi) D v - /q (z 2 ) + D K - /q (zi) D„+ /q (z 2 ) 



(Pu)i 



(21) 
(22) 

(23) 
(24) 

(25) 



5 



Pt. 



A N D n+/gf ( Zl ) A N D r+/9 r(z 2 ) + A N D T -/,t(«i) A N D n - /?f (z 2 )] _ (P L ) 



A' 



Pc = 



(Pi 



UN 



- (Pl)n 



Ar-/ g (zi) Ar-/g(z2)] 



(ft)i 



+ (ft) 



A UL,C, \ 1 

^12 (2l>22) = 



sm 



A UL,C, \ 



1 



cos z 

Zl ^2 



(ft - ft,c) 



(sin 2 2 ) 



7T Z 2 + Z 2 (1 + COS 2 6> 2 ) 



(ft - ft,c) ■ 



(26) 
(27) 
(28) 
(29) 



For fitting purposes, it is convenient to introduce favoured and disfavoured fragmentation functions, assuming 
in Eq. (TlOt): 



A N D 



P > TT+/u,d( Z ) 

Ar+/dt, g t(z,Rl_) 
P>Tr+/d,u(z) 



P>7T-/d,u( Z ) 



= 2M f l v (z)HP±) 



-p1/(p1) 



ApD 



AN P>7r-/uTJt(z,PA_) _ A N D 7r ± /s ^ s r(z, P± ) 



P > 7T-/u,d( Z ) 



ft 



7* 



2A^ s (z) h( PA 



-p1/(p!) 



ir(p 



(30) 
(31) 



with the corresponding relations for the integrated Collins functions, Eq. (17), and with A/" f ^ v dis (z) as given in 
Eq. jl2( with N q ° = Ng v>dis . 

We can now perform a best fit of the data from HERMES and COMPASS on A s ^ h+,ps) 
from the Belle Collaboration, on A^' C and An L ' C 



and of the data, 



Their expressions, Eqs. (14 1 and ( 25 )— ( 31 1 , contain the 
transversity and the Collins functions, parameterised as in Eqs. (|9])-(13). They depend on the free parameters 
a, (3, 7, 8, Nq , Nq and M^. Following Ref. [B] we assume the exponents a, j3 and the mass scale to be flavour 
independent and consider the transversity distributions only for u and d quarks (with the two free parameters 
Nu and iVJ). The favoured and disfavoured Collins functions are fixed, in addition to the flavour independent 
exponents 7 and 6, by iV f ^ v and . This makes a total of 9 parameters, to be fixed with a best fit procedure. 
Notice that while in the present analysis we can safely neglect any flavour dependence of the parameter f3 (which 
is anyway hardly constrained by the SIDIS data), this issue could play a significant role in other studies, like 
those discussed in Ref. l23l. 



II. BEST FITS, RESULTS AND PARAMETERISATIONS 

A. Standard parameterisation 

We start by repeating the same fitting procedure as in Refs. 6, 7J, using the same "standard" parameterisation, 
Eqs. (|6])-(13), with the difference that now we include all the most recent SIDIS data from COMPASS [TU] and 



HERMES [11] Collaborations, and the corrected Belle data [12] on A U 2 L and Af 2 c . Notice, in particular, that 
the A U 2 C data are included in our fits for the first time here. In fact, a previous inconsistency between A^ 2 L and 
A\ 2 data, present in the first Belle results [5T], has been removed in Ref. [T2"] . 

The results we obtain are remarkably good, with a total x 2 f of 0.80, as reported in the first line of Table |lj 
and the values of the resulting parameters, given in Table [H] are consistent with those found in our previous 
extractions. Our best fits are shown in Fig. [l] (upper plots), for the Belle A 12 data, in Fig. [2] for the SIDIS 
COMPASS data and in Fig. [3] for the HERMES results. 

We have not inserted the A§ Belle data in our global analysis as they are strongly correlated with the A\ 2 
results, being a different analysis of the same experimental events. However, using the extracted parameters 
we can compute the Aq L and Aq C azimuthal asymmetries, in good qualitative agreement with the Belle 
measurements, although the corresponding \ 2 values are rather large, as shown in Table [i] These results are 
presented in Fig. [l] (lower plots) . 

The shaded uncertainty bands are computed according to the procedure explained in the Appendix of Ref. [24) . 
We have allowed the set of best fit parameters to vary in such a way that the corresponding new curves have a 
total x 2 which differs from the best fit \ 2 by less than a certain amount A\ 2 - All these (1500) new curves lie 
inside the shaded area. The chosen value of A\ 2 = 17.21 is such that the probability to find the "true" result 
inside the shaded band is 95.45%. 
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TABLE I: Summary of the \ 2 values obtained in our fits. The columns, from left to right give the \ 2 P er degree of 
freedom, the total x 2 > an d the separate contributions to the total \ 2 of the data from SIDIS, AY2 , AY£ ', Aq L and Aq C . 
"NO FIT" means that the \ 2 for that set of data does not refer to a best fit, but to the computation of the corresponding 
quantity using the best fit parameters fixed by the other data. The four lines show the results for the two choices of 
parameterisation of the z dependence of the Collins functions (standard and polynomial) and for the two independent 
sets of data fitted (SIDIS, AY 2 L , AYf and SIDIS, Aq L , A% c ). 
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We have also performed a global fit based on the SIDIS and A a Belle data, and then computed the A12 
values. We do not show the best fit plots, which are not very informative, but the quality of the results can be 
judged from the second line of Table [ij which shows that although this time Aq L and Aq C are actually fitted, 
their corresponding \ 2 values remain large. This has induced us to explore a different functional shape for the 
parameterisation of Afq (z), Eq. ( [12] ), which will be discussed in the next Subsection. 

The difference between A12 and Aq is a delicate issue, that deserves some further comments. On the experi- 
mental side, the hadronic-plane method used for the extraction of Aq implies a simple analysis of the raw data, 
as it requires the sole reconstruction of the tracks of the two detected hadrons; therefore it leads to very clean 
data points, with remarkably small error bars. On the contrary, the thrust-axis method is much more involved as 
it requires the reconstruction of the original direction of the q and q which fragment into the observed hadrons; 
this makes the measurement of the A12 asymmetry experimentally more challenging, and leads to data points 
whith larger uncertainties. 

On the theoretical side, the situation is just the opposite: as the thrust-axis method assumes a perfect 



TABLE II: Best values of the 9 free parameters fixing the it and d quark transversity distribution functions and the 
favoured and disfavoured Collins fragmentation functions, as obtained by fitting simultaneously SIDIS data on the Collins 
asymmetry and Belle data on AY} and AY2 ■ The transversity distributions are parameterised according to Eqs. |9|, 
(111 and the Collins fragmentation functions according to the standard parameterisation, Eqs. (10 1, 1 12 1 and (131. We 
obtain a total xVd.o.f. = 0.80. The statistical errors quoted for each parameter correspond to the shaded uncertainty 
areas in Figs. [l]j3j as explained in the text and in the Appendix of Ref. [24] , 



N£ = 0.46t°-?° NJ = -1.00t£- 

g = 1-H±g:lg P = 3.64lj j° 

Ng v = 0.49t°; 2 ° iV£ s = -1.00t° 

7 = 1.06±°:« S = 0.07t ;f 7 
Ml = L50t 2 ;?2 GeV 2 
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FIG. 1: The experimental data on Ayt , A^-f (upper plots) and Aq L and Aq (lower plots), as measured by the Belle 
Collaboration [12] in unpolarized e + e~ — > hxh^X processes, are compared to the curves obtained from our global fit. 
The solid lines correspond to the parameters given in Table [II] obtained by fitting the SIDIS and the A\i asymmetries 
with the standard parameterisation; the shaded areas correspond to the statistical uncertainty on the parameters, as 
explained in the text and in Ref. [24] , Notice that the Aq L and Aq C data are not included in the fit and our curves, 
with the corresponding uncertainties, are simply computed using the parameters of Table [IT] 



knowledge of the q and q directions, the asymmetry can be reconstructed by a straightforward integration 
over the two intrinsic transverse momenta p±i and pj_2, transforming the convolution of two Collins functions 
into the much simpler product of two Collins moments [5], Eqs. (17) and (18). Instead, the phenomenological 
partonic expression of Aq requires more assumptions and approximations and is, on the theoretical side, less 
clean. 

One should also add that most of the large x 2 values found when computing A from the parameters of a 
best fit involving SIDIS and A12 data (or vice-versa) originate from the experimental points at large values of 
z\ or Z2 or both (see, for example the last points on the left lower panel in Fig. [IJ. Large values of z bring us 
near the exclusive process limit, where our factorized inclusive approach cannot hold anymore. 



B. Polynomial parameterisation 



In an attempt to fit equally well A\i and Aq (keeping in mind, however, the comments at the end of the 
previous Subsection) we have explored a possible new parameterisation of the z dependence of the Collins 
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FIG. 2: The experimental data on the SIDIS azimuthal moment A^j^ h+ ^ s ' as measured by the COMPASS Collabora- 
tion [TO] on proton (upper plots) and deuteron (lower plots) targets, are compared to the curves obtained from our global 
fit. The solid lines correspond to the parameters given in Table |TT| obtained by fitting the SIDIS and the A\2 asymme- 
tries with standard parameterisation; the shaded areas correspond to the statistical uncertainty on the parameters, as 
explained in the text and in Ref. |24j . 



function. Wc notice that data on A (z) seem to favour an increase at large z values, rather then a decrease, 



which is implicitly forced by a behaviour of the kind given in Eqs. ( 10 ) and ( 12 ) (at least with positive S values). 

In addition, an increasing trend of Aq(z) and An(z) seems to be confirmed by very interesting preliminary 
results of the BABAR Collaboration, which have performed an independent new analysis of e + e~ — > h± hi X 
data [25], analogous to that of Belle. 

This suggests that a different parameterisation of the z dependence of favoured and disfavoured Collins 
functions could turn out to be more convenient. Then, we try an alternative polynomial parameterisation 
which allows more flexibility on the behaviour of Afq (z) at large z: 



Nq (z) =N§ z[{l-a-b)+az + bz 2 



(32) 
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FIG. 3: The experimental data on the SIDIS azimuthal moment + as measured by the HERMES Collabora- 

tion are compared to the curves obtained from our global fit. The solid lines correspond to the parameters given 
in Table [IT] obtained by fitting the SIDIS and the A12 asymmetries with standard parameterisation; the shaded areas 
correspond to the statistical uncertainty on the parameters, as explained in the text and in Ref. [24] . 



with the subfix q = fav,dis, and — 1 < Nq < 1; a and b are flavour independent so that the total number of 
parameters for the Collins functions (in addition to Mh) remains 4. Such a choice fixes the term Nq (z) to 
be equal to at z = and not larger than 1 at z = 1. Notice that we do not automatically impose, as in 
Eq. ( fl2"| ), the condition |A/" g c (2;)| < 1; however, we have explicitly checked that the best fit results and all the 
sets of parameters corresponding to curves inside the shaded uncertainty bands satisfy that condition. 

We have repeated the same fitting procedure as performed with the standard parameterisation. When fitting 
the combined SIDIS, and Belle data, the resulting best fits (not shown) hardly exhibit any difference 
with respect to those obtained with the standard parameterisation (Fig. [l]). This can be seen also from the \ 2 ' s 
in Table [TJ where the third line is very similar to the first one. As a further confirmation, the corresponding 
best fit plots for A/" f ^ v dis (z), in case of the standard and polynomial parameterisations, plotted in Fig. |ij (left 
panel) practically coincide up to values of z very close to 1. 

The situation is different when best fitting the SIDIS data together with Aq L and Aq C ; in such a case 
the polynomial parameterisation allows a much better best fit, as shown in Fig. [51 upper plots. A reasonable 
agreement can also be achieved between the data and the computed values of Ayp and A\^ , as shown by the 
X 2 values in Table |lj and by the lower plots in Figj5] In this case the polynomial form of A/" f ° v dis (z) differs from 
the standard one, as shown in the right plots in Fig. [4] 

Notice, again, that the large \ 2 values of the computed A^2 is almost completely due to the last z bins, 
which correspond to the quasi exclusive region. Also, the larger \ 2 values corresponding to SIDIS data are 
mainly due to a slightly worse description of HERMES 7r~ azimuthal moments. The values of the parameters 
obtained using the polynomial shape of A/" f ° v dis (-z), Eq. (32), are given in Table 



III 



C. The extracted transversity and Collins functions; predictions and final comments 

Our newly extracted transversity and Collins functions are shown in Figs. [6] and [7J to be precise, in the left 
panels we show x Aj<q(x) — xhi q (x), for u and d quarks, while in the right panels we plot: 

zA N D h/qt (z) = zf d 2 p ± A N D h/q ,(z,p ± ) =zf d 2 P± ^ Hf*(z, P± ) = Az H^ 1/2)( '(z) (33) 
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FIG. 4: Plots of the functions (z) and JS fSJz ) for the favoured and disfavoured Collins functions as obtained by using 
the standard, Eq. ( 12 1, and polynomial, Eq. ( |32[ ) , parameterisations. On the left panel we show the results obtained by 
fitting the SIDIS data together with the A12 Belle asymmetries (both with standard and polynomial parameterisation), 
while on the right panel we show the corresponding results obtained by fitting the SIDIS data together with the Aq Belle 
asymmetries. 



TABLE III: Best values of the 9 free parameters fixing the u and d quark transversity distribution functions and the 
favoured and disfavoured Collins fragmentation functions, as obtained by fitting simultaneously SIDIS data on the Collins 
asymmetry and Belle data on Aq L and Aq° . The transversity distributions are parameterised according to Eqs. 
(111 and the Collins fragmentation functions according to the polynomial parameterisation, Eqs. ( 10 1, (32 1 and ( 13 1. We 
obtain a total xVd.o.f. = 1.01. The statistical errors quoted for each parameter correspond to the shaded uncertainty 
areas in Fig. |5j as explained in the text and in the Appendix of Ref. [24] . 



N£ = 0.36±g;il 
a = 1.06j 



H0.87 
-0.56 



i.oo.; n 

P = 3.66l 2 ..„ 



NJ = 



J 0.00 

3+5.87 



N c _ 1 nn+ 00 

cl — 2.36 q gg 

Ml = Omto'sj GeV 2 



N c — 

Jv dis 



-1 00+°' 19 
i.uu_ 00 



9 1 9+O.6I 
.12 



for h = tt^~ and q = u. The Collins results for d quark are not shown explicitly, but could be obtained from 
Tables [U and EED 



Fig. 



with the 



shows the results which best fit the COMPASS and HERMES SIDIS data on A s j^ h+0s) , together 



Belle results on and A^-f , using the standard parameterisation. The red solid lines correspond to 
the parameters given in Table [TTJ The shaded bands show the uncertainty region, which is the region spanned 
by the 1500 different sets of parameters fixed according to the procedure explained above and in the Appendix 
of Ref. [24]. The blue dashed lines show, for comparison, our previous results [7]: the difference between the 
solid red and dashed blue lines is only due to the updated SIDIS and AY% data used here, with the addition of 
, while keeping the same parameterisation. The present and previous results agree within the uncertainty 
band: one could at most notice a slight decrease of the new u quark transversity distribution at large x values. 

Fig.Rshows the results which best fit the COMPASS and HERMES SIDIS data on A s ^ h+<?is) , together with 
the Belle results on Aq L and A^ c , using the polynomial parameterisation. The red solid lines correspond to 
the parameters given in Table III This is not a simple updating of our previous 2008 fit [7] , as we use different 
sets of data (SIDIS and Aq rather than SIDIS and A12) with a different polynomial parameterisation. In this 
case the comparison with the 2008 results is less significant. If comparing the results of Fig.[6]and[7j one notices 
a sizeable difference in the favoured (u/ir + ) Collins function, and less evident differences in the transversity 
distributions. 

In Fig. H we show, for comparison with similar results present ed in Ref. , the tensor charge, corresponding to 
our best fit transversity distributions, as given in Tables ||l] and III Our extracted values are shown at Q 2 — 0.8 
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FIG. 5: The experimental data on Ag , Aq (upper plots) and A12 and A^ 2 (lower plots), as measured by the Belle 
Collaboration [12] in unpolarized e + e~ — > hxh^X processes, are compared to the curves obtained from our global fit. 
The solid lines correspond to the parameters given in Table [Till obtained by fitting the SIDIS and the Ao asymmetries 
with polynomial parameterisation; the shaded areas correspond to the statistical uncertainty on the parameters, as 
explained in the text and in Ref. [24] , Notice that the A^z and A^ 2 C data are not included in the fit and our curves, 
with the corresponding uncertainties, are simply computed using the parameters of Table |III| 



GeV 2 and compared with several model computations. One should keep in mind that our estimates are based 
on the assumption of a negligible contribution from sea quarks and on a set of data which still cover a limited 
range of x values. 

All other results are shown at the scale Q 2 = 2.41 GeV 2 . The evolution to the chosen value has been 
obtained by evolving at LO the collinear part of the factorized distribution and fragmentation functions. The 
TMD evolution, which might affect the k± and p± dependence, is not yet known for the Collins function. 
Consistently, it has not been taken into account for the other distribution and fragmentation functions. 

As BABAR data on An and Ao should be available soon, we show in Figs . |9| and [l0| our expectations, based 
on our extracted Collins functions. Fig. [9] shows the expected values of A^ 2 L , A^ 2 C , Aq and Aq C , as a function 
of Z2 for different bins of z±, using the parameters of Table \H\ obtained by fitting the SIDIS and the A12 Belle 



data with the standard parameterisation. Fig. 10 shows the same quantities using the parameters of Table |inj 
obtained by fitting the SIDIS and the Ao Belle data with the polynomial parameterisation. 

The Belle (and BABAR) e + e~ results on the azimuthal correlations of hadrons produced in opposite jets, 
together with the SIDIS data on the azimuthal asymmetry A^^ h+ ^ s \ measured by both the HERMES and 
COMPASS Collaborations, definitely establish the importance of the Collins effect in the fragmentation of a 
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FIG. 6: In the left panel we plot (solid red lines) the transversity distribution functions x hi q (x) = x Arg(i) for q 
with their uncertainty bands (shaded areas), obtained from our best fit of SIDIS data on A^ h+M and e+e 
on A12, adopting the standard parameterisation (Table |lTp . Simil arly , in the right panel we plot the corresponding first 
moment of the favoured and disfavoured Collins functions, Eq. (33 1. All results are given at Q 



dashed blue lines show the same quantities as obtained in Ref. [7] using the data then available on A 



= 2.41 GeV^. The 

•taWh+fcO and a ul 
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transversely polarised quark. In addition, the SIDIS asymmetry can only be observed if coupled to a non negligi- 
ble quark transversity distribution. The first original extraction of the transversity distribution and the Collins 
fragmentation functions [6j[7], has been confirmed here, with new data and a possible new functional shape of 
the Collins functions. The results on the transversity distribution have also been confirmed independently in 
Ref. 0. 

A further improvement in the QCD analysis of the experimental data, towards a more complete understanding 
of the Collins and transversity distributions, and their possible role in other processes, would require taking into 
account the TMD-evolution of Arq(x, k±) and A D h i q -\ (z,p±). Great progress has been recently achieved in the 
study of the TMD-evolution of the unpolarized and Sivers transverse momentum dependent distributions 33 37 
and a similar progress is expected soon for the Collins function and the transversity TMD distribution [38 . 
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FIG. 9: Estimates, obtained from our global fit, for the azimuthal correlations A" 2 , A{ 2 , Aq an d Aq in unpolarized 
e + e~ — > h\h 2 X processes at BaBar [55]. The solid lines correspond to the parameters given in Table [TTj obtained by 
fitting the A\i Belle asymmetry; the shaded area corresponds to the uncertainty on these parameters, as explained in 
the text. 
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FIG. 10: Estimates, obtained from our global fit, for the azimuthal correlations A^^ , A\ 2 , Aq L and Ag C in unpolarized 
e + e~ — > h\h 2 X processes at BaBar [25]. The solid lines correspond to the parameters given in Table III obtained by 



fitting the Ao Belle asymmetry; the shaded area corresponds to the uncertainty on these parameters, as explained in the 
text. 



